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Nematic textures in spherical shells. 
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The equilibrium texture of nematic shells is studied as a function of their thickness. For ultrathin 
shells the ground state has four short ^ disclination lines but, as the thickness of the film increases, a 
three dimensional escaped configuration composed of two pairs of half-hedgehogs becomes energeti- 
cally favorable. We derive an exact solution for the nematic ground state in the one Frank constant 
approximation and study the stability of the corresponding texture against thermal fluctuations. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

The study of liquid crystal phases benefits from geo- 
metrical reasoning in two important ways. Firstly, liquid 
crystal elasticity can often be cast in terms of the cur- 
vature of equipotential lines (or surfaces) that map out 
the corresponding textures. Second, the observed tex- 
tures are strongly affected by geometric and topological 
constraints imposed by the presence of boundaries con- 
fining the system. The liquid crystal ground state results 
from the competition between the energetic requirement 
of minimizing the " curvature of the texture" and the geo- 
metric frustration introduced by boundaries that impart 
a preferred curvature at the edge of the sainple that often 
cannot propagate across the system [ll S 13 ■ 

The boundary conditions can be controlled experimen- 
tally with the possibility of designing molecular systems 
with intriguing technological applications 3 • For exam- 
ple, colloidal particles coated by a very thin nematic layer 
have in their ground state four disclinations sitting at the 
vertices of a tetrahedron. Each coated colloidal particle 
can then be viewed as the fundamental building block of a 
self assembled lattice with tetravalent coordination. The 
" bonds" between the colloidal particles could be provided 
by chemical linkers attached at the four "bald spots" at 
the cores of the four disclinations present in each colloid 
13 • A second example, is provided by self-assembled sys- 
tems of block copolymers Q which are a promising tool 
for "soft lithography" on both flat and curved substrates 
Q. In addition, liquid crystals in confined geometries 
provide an arena for physicists and mathematicians in- 
terested in applications of geometrical and topological 
ideas to material science ^ i^ |1^ . 

In this work we present a theoretical study of liquid 
crystal phases (focusing on vector, nematic and hexatic 
order) confined in a spherical shell of varying thickness 
with the director assumed to be tangent to the two in- 
terfaces. We first consider the two dimensional regime 
where a nematic film coats a quenched spherical surface 
such as a colloidal particle in solution or the interface of, 
say, a water droplet in oil. The presence of topological 
defects in the ground state for ordered states on spherical 
surfaces is unavoidable [ll|, [l^, [l3 ■ Recent experimen- 
tal and theoretical investigations of spherical crystallog- 



raphy have provided an alternative context to study the 
constraints posed by the compactness of the underlying 
curved space [ij,[ia|- More recent explorations have con- 
centrated on 2D ordered phases confined to interfaces of 
varying Gaussian curvatur e Il6lll7lll8| | as well as dynam- 
ically fluctuating surfaces [I^ |23 . 

As the thickness of a nematic film increases, an es- 
caped three dimensional texture, also strongly influenced 
by the spherical topology and the boundary conditions, 
become energetically favored with respect to planar tex- 
tures. This instability destabilizes the tetravalent ne- 
matic texture on colloids. In this paper, we estimate the 
thickness of the nematic film above which a texture with 
four radial disclination lines of charge s = ^ becomes un- 
stable to four half-hedgehogs. The two competing tex- 
tures studied in this paper are shown in Fig. ^ We 
also discuss the possibility of hysteresis between the two 
textures. 

The organization of this paper is as follows. In sec- 
tion ^] we derive exact solutions for the ground state of 
spherical films of tilted molecules and nematogens within 
isotropic elasticity by using the method of conformal 
mappings. In the notation of References 0, [i3j these 
situations correspond to order parameters described by 
a bond angle with p = 1, 2 fold symmetry in the tangent 
plane of the sphere (see Appendix El- A mathematical 
justification for our approach is provided in Appendix IbI 
where the same technique is illustrated in the context of a 
more familiar flat space problem. In section UTTl we study 
the stability of liquid crystal textures to thermal fluctua- 
tions by means of a normal mode analysis whose details 
are relegated to Appendices El a-nd [0 The stability of 
the valence-four texture against escaped solutions is con- 
sidered in section IIVI where a phase diagram is derived 
with the thickness as a control parameter. The texture 
distortions caused by the elastic anisotropy between bend 
and splay deformations are briefly considered in section 

m 



II. TEXTURES 

The liquid crystal free energy for molecules embedded 
in an arbitrary frozen surface can be written in the one 








FIG. 1: (Top panel) Two-dimensional texture characterized 
by four short disclination lines at the vertices of a tetrahe- 
dron inscribed in the sphere. The surface texture shown (in- 
scribed on the surface of a baseball) is invariant throughout 
the thickness of the shell. (Bottom panel) Cut view of the 
escaped three dimensional texture given by two pairs of half 
hedgehogs located at the north and south poles of the sphere. 



constant approximation as 



F= — f dAD^n^{u)D'nj{u) , 



(1) 



where u = {^1,1*2} is a set of internal coordinates, n(u) 
is the liquid crystal director defined in the tangent plane, 
Di is the covariant derivative with respect to the metric 
of the surface and dA is the infinitesimal surface area 
[S IM El HI. The free energy of Eq.© is invariant 
upon rotating each molecule n(u) by the same (arbi- 
trary) angle with respect to any axis of rotation per- 
pendicular to the local tangent plane. The treatment of 
systems with a p-fold symmetry is straightforward pro- 
vided that the one Frank constant approximation is used 
for p = 1 and p = 2 and the consequences of any addi- 
tional couplings to curvature neglected [23|. This choice 
of free energy implies that the minimal energy configu- 
ration will be given locally by neighboring n(u) vectors 
which differ only by parallel transport. The curvature 
of the surface induces "frustration" in the texture. In 
fact, by Gauss' "Theorema egregium" [g, tangent vec- 
tors parallel transported along a closed loop are rotated 



by an amount equal to the Gaussian curvature integrated 
over the enclosed area. On a sphere, this theorem in- 
sures that the nemat ic g round state always has four ex- 
cess disclinations [ll|, [13 . More generally, the sum of the 
topological charges on any closed surface is equal to the 
integrated Gaussian curvature, implying a minimum of 2 
and 6 disclinations in the ground state of tilted molecules 
and hexatics, respectively. 

We introduce a local angle field a{u), corresponding to 
the angle between n(u) and an arbitrary local reference 
frame, we can rewrite the free energy introduced in Eq.lQJ 
as: 



F 



1 



K / dSg'^{d,a - Ai)idja - A^) 



(2) 



where dS = d^u^/g, g is the determinant of the metric 
tensor gij and Ai is the spin-connection whose curl is the 
Gaussian curvature G'(u) H, |23. On a sphere of radius 
R parametrized by polar coordinates (0, 0), the only non 
vanishing components of the (inverse) metric tensor are 



9' 



B.sin0 



and 5* 






A convenient choice of of 



the spin connection (which plays the role of the vector 
potential) is discussed in Appendix ^ The simplified 
free energy in Eq.(|2J) is the starting point of our analysis. 



A. Tilted molecules on a sphere 

The orientational order of molecules tilted by a con- 
stant angle with respect to a spherical interface can be 
modelled by a vector field n{6, (j)) defined in the local 
tangent plane on which the molecule has a fixed length 
projection [l^. To determine the ground state of the liq- 
uid crystal texture, we minimize the Frank free energy of 
Eq.®. As discussed above, the topological charges must 
sum up to 47r, the integrated Gaussian curvature of the 
sphere |al23- For a vector field {p = \) the texture with 
only two defects of charges -|-27r minimizes the Frank free 
energy and satisfies the topological constraint. Since the 
defects repel each other they preferentially sit at two an- 
tipodal points that we can designate as the north and 
south pole of the sphere. If the splay and bend cou- 
pling constants of the nematic are equal, then there is a 
large degeneracy in the ground state arising from the in- 
variance of the vector free energy in Eq. Q under global 
rotations a(u)^a(u)-|-c, where u= 9,ip. One represen- 
tative texture is a "sink" and a "source" of n(u) at the 
two poles. In this splay rich texture n(u) is parallel to 
the lines of longitude on a sphere. In a bend rich tex- 
ture, related to the previous by a -I rotation about the 
local normal to the surface, n(u) is everywhere parallel 
to the lines of latitude. Any other rotation of n(u) that 
makes an arbitrary constant angle with respect to this 
texture is an acceptable solution for the ground state of 
the molecules. 

As we now show, this degeneracy is lifted when K3 ^ 
Ki. Indeed the effect of distinct splay and bend elastic 
constants Ki and Kj, (the twist elastic constant K2 is 



absent in two dimensions) is to select the bend-rich tex- 
ture if Ki > K^ or the splay-rich one if K^ > Ki . The 
intermediate configurations obtained by a global rotation 
of the director are now unstable. Assume for simplicity 
that K3 > Ki. In this case, it is convenient to recast the 
Frank free energy (see Appendix A) as follows: 

+ {K3 - K,){T) X n)^] , (3) 

where the covariant derivatives is expressed in terms of 
the Christoffel connection, F^^ , 



F^ 

-^ it 



(4) 



(5) 



Di n' ^ di n^ 
and the covariant form of the curl squared is |2S 

(V X n)2 = {D, Uj - Dj n,) {D' n^ ~ D^ n' 

The first term in Eq.Q resembles the Frank free en- 
ergy in the one coupling constant approximation and is 
minimized by choosing the sink-source (or "lines of lon- 
gitude" ) solution. The second term (which is positive 
definite) will vanish for this texture since the sink-source 
texture is bend free. All other textures have a higher 
energy. 

A similar argument can be used to prove that the two 
vortex-configuration which follows the lines of latitude 
is the minimum of the free energy when Ki > K^ by 
rewriting the Frank free energy as 

F=]^j(f^^ [A'3 (A n' ) ( D' n,) 

+ (K, - K3) (D • n)2] , 

where the covariant form of the divergence reads 

1 



(6) 



D 



V9 



dr {^n^) 



(7) 



The latitudinal texture minimizes the first term of Eq. @ 
while the second vanishes because this texture is splay 
free. Any deviation from the splay-free latitudinal tex- 
ture will only increase the energy. 

The energy of both textures can be expressed as a func- 
tion of the anisotropy parameter, e, and the mean of the 
elastic constants, AT, 



K = 



K3 






(8) 
(9) 



and the radius of the sphere, R, scaled by the short dis- 
tance cutoff a. The resulting free energy for arbitrary e 
reads (see Ea. (|A16|l ') 



F = 27rA:(l- |e |) In 



0.3 



(10) 



The conclusions of this section are summarized in Fig. 
El which suggests that there is a discontinuous first or- 
der transition when e passes through zero. This analysis 
mirrors similar arguments valid in the plane |25j. 





-1 



FIG. 2: Schematic illustration of the phase diagram of the 
texture of tilted molecules on a sphere as a function of the 
anisotropy parameter e superimposed on a plot of the free 
energy, F, stored in the texture versus e. The two competing 
ground states (top view) are characterized by either pure bend 
(lines of longitude configuration at left) or pure splay (lines 
of latitude configuration at right). 



B. Nematic Texture 

The nematic texture of very thin spherical shells of 
nematic liquid crystal with tangential boundary condi- 
tions can be analyzed within the one Frank constant 
approximation by using the method of conformal map- 
pings whose mathematical justification is illustrated in 
Appendix IbI bv means of a simpler example. 

An elegant argument introduced by Lubensky and 
Frost in Ref. 13] shows that the ground state of nemato- 
gens on a sphere is given by 4 disclinations of topological 
charge s = 1/2 sitting at the vertexes of a tetrahedron. 
The energy of single disclinations is proportional to the 
square of its strength. As a result, the longitudinal and 
latitudinal textures derived for tilted molecules in section 
III Al are unstable since their energies can be lowered by 
splitting each ,s = 1 defect at the north and south pole 
into two s = 1/2 disclinations and letting them relax to 
their equilibrium positions at the vertexes of a tetrahe- 
dron where they are as far away from each other as pos- 
sible. According to a calculation in Ref. [13, the energy 
Fg of a sphere of radius R with in plane orientational or- 
der and 2n interacting minimal disclinations for a p-fold 
order parameter is given by 



R = 2ttK h 



p \ a 



(11) 



where the {cp} are constants depending on the symmetry 
of the order parameter and the defect core energy while h 
is the thickness of the liquid crystal layer. The numerical 
values of the relevant constants are ci = and C2 — —0.2. 
When p = 2 is chosen in Ea. (|ll|) the elastic energy is 
indeed smaller than the corresponding value for p = 1 in 




FIG. 3: Schematic illustration of the baseball texture of a 
thin nematic shell. The same texture is reproduced from a 
different perspective in the top panel of Fig. 




FIG. 4: Graphical construction of the stereographic projec- 
tion. Regions close to the north pole have larger images in the 
conformal plane than regions of equal areas close to the south 
pole. The stereographic projection preserves the topology of 
the surface provided all points at infinity are identified with 
the north pole. 



the limit R ^ a, in agreement with related arguments 
given in Ref.|5l|. 

To obtain an algebraic expression for the texture we 
proceed as illustrated in Appendix El and seek a func- 
tion fl{x, y, z) = ^{x,y, z)+i"i'{x, y, z) which is harmonic 
on the sphere except for two arcs connecting the defects 
in pairs. The calculation for nematogens described be- 
low was suggested to us by F. Dyson [23. The function 
$(x,2/, z), which we can interpret as an electrostatic po- 
tential, takes equal and opposite values on the two arcs 
and is equal to zero on a baseball-like seam (see Fig. 
^ which divides the sphere into two congruent regions. 
The nematic director is then oriented (up to a global ro- 
tation) along the contour lines of <I>(a;, y, z), that is the 
equipotential lines of this " curved space capacitor" . In 
this analogy, the contour Hues of '^{x,y,z) are electric 
field lines, hence they correspond to a valid texture where 
the director is rotated locally by ^ with respect to the 
equipotential lines. The arcs can be either great-circle 
arcs extending more than half way round the sphere or 
short great-circles arcs connecting the same pair of de- 
fects along the shortest path. The first choice leads to 
equipotential lines whose seam resembles in shape that 
of a baseball. If the second choice is made the pattern of 
equipotential lines would not deviate much from concen- 
tric circles and the seam would look more like the seam 
of a cricket ball. We will explicitly show that the two 
choices are equivalent since the equipotential lines of the 
first solution are field lines of the second and vice versa. 

We choose the arcs connecting the defect pairs along 
great circles and we take the four defects labelled by A, 
B, C, D to lie at the vertices of a tetrahedron inscribed 
on a sphere of radius 1 and whose north and south poles 
are N = (0, 0, 1) and S = (0, 0, -1) respectively 






We now perform a stereographic projection (see Fig. 0)) 
that maps every point on a unit sphere centered on the 
origin onto the plane z = — 1 according to the rule 



(13) 



The coordinates of the image points (connected to points 
on the sphere by dashed lines in Fig. 0)) are given by 

2x 




b = 



1- z 
2?/ 



1 



(14) 



Upon transforming to a complex coordinate w = a + ih, 
the four tetrahedral points of Eq. \i2\ are mapped onto 

^' = p{l + i) , B' =p{-l-i) , 

C = q{-\ + i) , D' = q{\-i) , (15) 



where 



P 



Vi + l 



q 



\/3-l 



(16) 



(In this section p does not refer to the symmetry of the 
order parameter). The great arc passing through the 
south pole (corresponding to one capacitor plate in the 
electrostatic analogy) maps onto the segment A'B', as 
illustrated schematically in the top panel of Fig. while 
the great arc through the north pole maps onto the two 
semi-infinite segments of the line a + b = which bracket 
CD' . We can fold the two cuts in the w plane back on 
top of each other by mapping the w plane onto the u 
plane via 



u = —iw 



(17) 



As shown in the middle panel of Fig. [S] the images A and 
B of A' and B' now both lie on the real axis at 2p^ while 
the images of C and D' now lie at —2q^. The two cuts in 
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FIG. 5: Illustration of the change in the branch cut structures 
of the complex function i}{v) describing the nematic texture 
after performing a series of conformal transformations. In 
the top panel we have simply performed a stereographic pro- 
jection from the sphere. The middle panel shows the "fold 
up" transformation of Eg. 1171 whereas the bottom panel cor- 
responds to the transformation in Ea. ll2H that symmetrizes 
the positions of the cuts. 



the u plane are both on the real axis, running from zero 
to 2p^ and from — 2g^ to minus infinity. On the sphere, 
these correspond to geodesies connecting defects which 
stretch more than halfway around the sphere. In order to 
make the cuts symmetric with respect to the imaginary 
axis (see the bottom panel of Fig. [SJ we search for a 
conformal transformation that maps the following four 
points in the complex w-plane to four points on the real 
axis of a complex w-plane 



Mo = - 


-^ uo = 


--k , 


Ml == -2q^ - 


-> Ul = 


= -k , 


M2 = -OO - 


-^ V2 = 


= -1 , 


M3=2p2 - 


-^ U3 = 


= 1 . 



(18) 



In order to fully determine the conformal transformation 
we need to determine the value of k. This can be done 
by using a standard relation in the theory of conformal 
transformations ^3 



Mo - Ml M3 - M2 Mo ~Vl M3 - V2 



Mo - M2 M3 



Ml 



Mo - M2 M3 



Ml 



(19) 



Upon inserting the points of Ea. (|18|l into Ea. (|19|) . we 
determine the value of k (less than one) 



k = 



2V2-P 

P + 2V2 



(20) 



Equation l|18|) contains four independent relations so we 
are still left with three conditions to determine the three 
independent coefhcients {a, /9, 6} of the bilinear confor- 
mal transformation that implements the mapping illus- 
trated pictorially in the bottom plate of Fig. El 



(21) 



au + P 



The required coefficients needed to implement the map- 
ping in Ea. H18|) are 



a = — 1 , 

/3 = 2p{2V2+p) 

5 = 2p (2\/2 - p) 



(22) 



To solve Laplace's equation, we desire a function Q.{v) 
which is analytic except on the two cuts on the real axis, 
and whose real part takes constant values on the cuts. By 
symmetry, rL{v) is an odd function of v, and its real part 
$(m) is zero when the real part of v is zero. Therefore the 
image of the seam in the v plane is simply the imaginary 
axis SRe m = 0. Upon substituting for v using Equations 
(|21|l and (|17|) . the condition 5Re m = becomes 



16 + VQto(m;'^)- |w|'' 







(23) 



With the help of Equations H14fl and (|13|l . we can now 
write down the equation of the seam explicitly in the 
original cartesian coordinates [2a| 



z^{2 + Vi)xy , 
or in spherical polar coordinates {0, 9} as 
cos 9 



sm 



1 + — ] sin 20 



(24) 



(25) 



The seam defined by the line of zero potential, is repre- 
sented for different orientations of the sphere in Fig. 
Its contour length Z, on a unit radius ball, is readily calcu- 
lated upon integrating the expression for the infinitesimal 
arc of the seam 



dl 



Uin^f^] +1 d9 



\d9 J 



(26) 



from 9min ~ 0.69 radians to 9„iax ~ 2.44 radians and 
multiplying the result by four in view of the symmetry 
of the seam. The values of 9min and 9max are obtained 
from Eq. H25|) by setting equal to ^ and ^ respectively. 
Upon using Ea. lf^ to substitute 0(6') in Eg. lf^ . we 
obtain I ~ 9.09 for a sphere of unit radius. The seam 






FIG. 6: DifTerent views of a track of parallel nematogens 
which partitions the sphere into two equal areas, each contain- 
ing two s — ^ disclination defects. It resembles a "fattened" 
version of the seam of a baseball. 



is longer than the equatorial circumference by slightly 
less than 50%. 

The branch cut structure in the v plane is sufficiently 
simple to allow a guess of the corresponding analytic 
function f2(ti). A function with cuts from A; to 1 and 
—k to —1, whose real part is equal and opposite on the 
two cuts and with a single imaginary period around any 
curve separating the cuts is easily identified to be a stan- 
dard elliptic integral 



n{v) ^ / [{k^ - t^){i - 1^)] ^dt 



(27) 



with V given in terms of w by Equations (fTTjl and (j?T|l . 
The nematic director is oriented (up to a global rotation) 
along the contour lines of the imaginary or (real part) 
of i^{v). The equipotential (red) and field lines (black) 
of fl(v) are conveniently plotted using the stereographic 
projection plane w — a + ib in Fig. [7\ along with the 
positions of the disclinations (green dots). It is easy to 
switch from the stereographic-projection plane w = a+ib 
of Fig. [3 to spherical polar coordinates {9, (f) by using 
the relation (reviewed in Appendix IE)) . 



w = 2Rcot \ - I e 



itf) 



(28) 



where R is the radius of the sphere. 

If we had constructed the base-ball with cuts along the 
short geodesies connecting the defects, then the form of 
the texture given in Ea. (|27|l would be the same, but the 
parameter k in the elliptic integral would be given by 

fc-(V3-V2)2(V2 + l)2 , (29) 

instead of Eq. H2UI) . The equation for the seam becomes 

z = -{2 ~ ^/i) xy (30) 




FIG. 7: Illustration of the nematic texture in stereographic 
projection (Color online). The red and black field lines corre- 
spond to two energetically degenerate (in the one Frank con- 
stant approximation) families of bend and splay rich textures. 



The dots indicate the four tetrahedral s 
fects. 



disclination de- 




and the corresponding equipotential and field lines are 
plotted in black and red respectively in Fig. The 
expression reads 



(31) 



in spherical polar coordinates, and leads to the same con- 
tour length I as Eq. (|25|l . Thus the two choices of arcs lead 
to two equivalent textures differing only by a ^ about 
the local normal. As in Section FlI Al the degeneracy in 
energy between the red and black flow lines in Fig. [3 
is lifted upon considering the effect of elastic anisotropy 
generated by a different energy cost for bend and splay. 



III. STABILITY OF LIQUID CRYSTAL 
TEXTURES TO THERMAL FLUCTUATIONS 

In this section, we study the stability of liquid crystal 
ground states to thermal fluctuations [3 . To explore the 
fidelity of directional bonds at finite temperatures, we 
employ a Coulomb gas representation of the liquid crystal 
free energy (in the one Frank constant approximation) 
obtained by substituting in Eq.jJl the relation 



I'^^do.idpO - Ap) = s(u) - G(u) = n(u) 



(32) 



where ^°'^ 

m, G(u) 



is the covariant antisymmetric tensor 

is the Gaussian curvature and s(u) = 

VLi) is the disclination density with A^^ 



defects of charge qi at positions u^ . The final resuh is an 
effective free energy whose basic degrees of freedom are 
the defects themselves [HlllEJ: 



F^§ldA 



fdA fdA' 7i(u) r(u,u') ?i(u') 



(33) 



The Green's function r(u, u') is calculated (see Appendix 
IX| by inverting the Laplacian defined on the sphere 



r(u,u') 



(34) 



and we have suppressed for now defect core energy con- 
tributions which reflect the physics at microscopic length 
scales. Equations (|32|) and H33|) can be understood by 
analogy to two dimensional electrostatics, with the Gaus- 
sian curvature G{u) (with sign reversed) playing the role 
of a uniform background charge distribution and the 
topological defects appearing as point-like sources with 
electrostatic charges equal to their topological charge qi. 
The charge q can be defined by the amount 9 increases 
along a counterclockwise path enclosing the defect's core. 
On a generic surface, the defects tend to position them- 
selves so that the Gaussian curvature is screened: typ- 
ically, the positive ones are attracted to peaks and val- 
leys while the negative ones to the saddles of the sur- 
face [13, [13. This geometric potential is ruled out by 
symmetry on an undeformed sphere since the Gaussian 
curvature is constant. The Gaussian curvature plays the 
role of a uniform background charge fixing the net charge 
of the defects consistent with the topological constraint 
imposed by the Poincare-Hopf theorem (see section Hll 
and References |Ml23)- The equilibrium positions of the 
defects are then determined only by defect-defect inter- 
actions which are proportional to the logarithm of their 
chordal distance (see Appendix ^ according to 



F ■ 



irK 



y UiTij In [1 — cos f3i. 



(35) 



ii^j 



where the integers Ui and rij describing the singularities 
associated with each defect, and the integer p controls 
the period — of the orientational order parameter. The 
"topological charge" describing the rotation of the order 
parameter around each defect is given by Sj = — . The 
geodesic angle (3ij subtended by the two defects at posi- 
tions Ui — {9i, (/)[} and Uj = {9j, 0j} can be conveniently 
recast in terms of their spherical polar coordinates 



C0S/3r 



COS f , cos ( 



sm f , sm ( 



(36) 



As a simple example, we first consider the case of a 
Z = 2 colloidal particle (p=l) with two antipodal defects 
of index 1. We can study the effect of thermal disruption 
of the ground state by setting (3ij = n + 9 in Ea. l|5S|l 
and expanding in the bending angle 9. The resulting free 
energy, apart from an additive constant, reads 



F: 



nK 



(37) 



Upon applying the equipartition theorem we obtain in 
the hmit K > fc^T 



(cos/3y) 



1 



-1 + - 
-1 



ksT 



ttK ' 



(38) 



which describes the fidelity of tt antipodal " bonds" of a 
divalent colloidal particle. 

The effect of thermal fluctuations on the tetrahedral 
ground state of nematic molecules confined on the sphere 
{p—2 and all Sj — 2) can be studied by means of a nor- 
mal mode analysis. The basic results sketched in were 
obtained by a slightly different method. Here we describe 
an alternative treatment in some detail and extend our 
analysis to hexatic and tetratic defect arrays (see Ap- 
pendix Q . 

We start by defining a generalized array of defect coor- 
dinates {qi} as a 27V dimensional vector, where N is the 
number of defects in the ground state (or, equivalently, 
the valence of the colloidal molecule) . iV = 4 in the case 
of the tetrahedron. The first A^ entries of the vector q 
are the longitudinal deviations of the N defects from a 
perfect tetrahedral configuration while the remaining N 
components describe defect displacements along the lines 
of latitude of a sphere of unit radius. As a result, the de- 
viations of the i defect from its equilibrium configura- 
tion {9^ , (/)/'} are parameterized by the two independent 
components of the vector q 



qN+i ~ 



'n- 



(39) 



The relations in Eq. H39() can be used to reexpress Eq. H36() 
in terms of the components of the displacements vector 
Qj, with the result. 



cos Pij — cos (61" -I- qi) cos (0/ + qj) 



sin {91' + q,) sin {91' + qj] 



cos 



qN+i 

sin 9P 



iN+j 



Upon substituting Eg. 1)401) in Ea. (|35() . the free energy F 
can be expanded around the equilibrium configuration 
to quadratic order in qi with the result (apart from an 
additive constant) 



F 



Je^^^ 



ij qi qj 



(41) 



where the matrix, M^ , describing the deformation of the 
tetrahedral molecule is naturally defined as 



My- 



d^F 
dqi dqj 



(42) 



qi,q j=0 



The eigenvalues of this matrix can be classified according 
to the irreducible representation of the symmetry group 



(40) 
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TABLE I: Character for the irreducible representations of the 
tetrahedral point group together with the character of the 
eight-dimensional representation E generated by the defect 
displacements of a tetravalent colloid. 



3d 


E 


8C3 


3C2 


654 


60-d 


Al 


1 


1 


1 


1 


1 


A2 


1 


1 


1 


-1 


-1 


E 


2 


-1 


2 








Fi 


3 





-1 


1 


-1 


F2 


3 





-1 


-1 


1 


J] 


8 


-1 












of the tetrahedron; their degeneracies can be determined 
purely from the group theoretical relation |29l l3(]| 



^(7) 



1 Y^ (^)* (S) 

-} ,9iXi Xi 



(43) 



where n^'^' is the number of frequency degenerate nor- 
mal modes that transform like the irreducible represen- 
tation labeled by 7, gi is the number of symmetry op- 
erations of the tetrahedral point group in the i*^ class, 
g = ^^ gi = 24 is the total number of symmetry opera- 
tion in the group, Xi is the character of the i*'* class in 

the irreducible representation labelled by 7 while xl is 
the corresponding character for the reducible representa- 
tion formed by the defects' displacements. 

The information necessary to apply Ea. (|43() to a 
tetravalent colloid is collected in Table |I| The top row 
contains the five symmetry class {E, C3, C2, iS'4, ad} con- 
tained in the tetrahedral point group 3^, correspond- 
ing respectively to the identity, three and two fold rota- 
tions, four fold rotatory-reflections and reflection through 
a plane of symmetry |29l | . The number of symmetry op- 
erations gi included in the i*'' class also appears in the 
top row: thus, {gi} = {1,8,3,6,6} where the same or- 
dering used above to list the classes has been adopted. 
The left most column of Table Q] lists the one, two and 
three dimensional irreducible representations of the tetra- 
hedral group {Ai,A2, E, Fi,F2}, along with the eight di- 
mensional representation E generated by the defect dis- 
placements. The entries of the table list the characters 
corresponding to each class of the five irreducible repre- 
sentations, Xj , and in the last row the corresponding 



,(7) 

-i ' 



characters, Xi : for the eight dimensional representa- 
tion. The former are tabulated from standard group the- 
oretical treatments while the latter needs to be worked 
out from the traces of the transformation matrices that 
describe how the displacement coordinates qi transform 
under the action of each symmetry element in the group. 
These manipulations are rather cumbersome, especially 
for the "icosahedral molecule" arising when a spherical 
surface is coated with a pure hexatic layer (see Appendix 

In the rich literature on molecular vibrations a set 



of empirical rules has been developed to write down 
the characters by examining only the transformation of 
the three dimensional cartesian displacements of the few 
atoms whose equilibrium positions are not altered by the 
symmetry operation. In Appendix |C1 we provide analo- 
gous rules that simplify the task of finding the x} char- 
acters by incorporating the constraint that each atom is 
confined on a sphere and hence only two orthogonal dis- 
placements need to be considered as shown in Ea. (|S^ . 

The interested reader is referred to Appendix [CI for 
a more comprehensive mathematical justification of the 
normal mode analysis applied to the tetrahedral colloid 
and to the more complicated cases of hexatic Z = 12 and 
tetratic order Z = 8. Here, we simply summarize the 
results of applying Ea. H43|) in conjunction with Table Q] 
to find the degeneracies of the eigenvalue spectrum of the 
matrix Mij. The representation S contains (only once) 
the three dimensional representations F2 and Fi as well 
as the two dimensional representation E. 



Y.^F2+Fi+E 



(44) 



The three normal coordinates with vanishing frequency 
correspond to the three rigid body rotations and belong 
to the Fi irreducible representation |29l |3Q| . We are left 
with a doublet (E) and a triplet (F2) corresponding re- 
spectively to two shear-like twisting deformations of the 
tetrahedron and to three stretching and bending modes 
of the cords joining neighboring defects. 

This symmetry analysis is confirmed by direct diago- 
nalization of the matrix Alij which leads the following set 
of eigenvalues Xi 



{AJ 



SnK 



{0,0,0,1,1,2,2,2} 



(45) 



In Sectional we also list the eigenvectors Wi oi Mij. The 
displacement coordinates are readily expressed in terms 
of the eigenvectors 



'1^ = U^'' 



'3 ' 



(46) 



where the unitary matrix U diagonalizes M and hence 
the free energy of Eq. (|41|l and is defined by 



U MU^ 



Diag {Xi 



(47) 



Its construction is easily achieved by the standard Gram- 
Schmidt orthogonalization procedure to the eigenvectors 
{wi} = {wi, ...,ws}, where the same ordering chosen in 
listing the eigenvalues in Ea. (|45|l is implicitly assumed. 
The resulting orthogonal basis vectors are the rows of the 
8x8 matrix U. 

We are now in a position to evaluate (cos /3y ) where the 
thermal average is performed with the Boltzman weight 
obtained from the free energy in Eq. (|41|l which is now di- 
agonal. Note that for the tetrahedron any choice of pair 
of defects labelled by i and j (where i ^ j) will lead to 
the same answer, unlike the less symmetric cases of the 
twisted cube (p=4) and the icosahedron (p=6) consid- 
ered in AppendixlO The bending angle cos (3ij in Eq. 14U|) 



can be Taylor expanded in the qi. The resulting expres- 
sion is rather cumbersome, but once the displacements 
{qi\ are reexpressed in terms of the normal coordinates 
{wi} (by means of Ea. (|l^ ). cos (iij reduces to 



COS^y = -■; 



1 



(48) 



where the only eigenmodes {'WQ,W'j^wg,} appearing in 
Ea. (|48|l correspond to the bending triplet of Ea. (|45f) . 

It is now easy to perform the thermal average by Gaus- 
sian integration of the energetically degenerate eigen- 
modes, with the result Q 



(cos/3i. 



(( 



w, 



6/ 



1 2 

1 IQkBT 
3^ Q-kK 



oD) 



(49) 



IV. VALENCE TRANSITIONS IN THICK 
NEMATIC SHELLS 



FIG. 8: Side view of the 2D nematic texture ansatz solution 
in Ea. l541 . The cylindrical slab has height h and radius R. 
The arrows can be interpreted as flow lines of a fluid entering 
a narrow and long channel from a point source located on the 
top plate. To obtain the nematic texture the vectors need to 
be normalized and viewed as rods with both ends identified 
as shown in Fig. |3] 



In this section, we study the crossover from a two to 
three dimensional regime as the thickness of the spherical 
shell, h, increases. For thicker shells, three dimensional 
defect configurations ("escaped" in the third dimension) 
compete with the planar textures described in the pre- 
vious sections, leading to a structural transition and a 
change in valence from Z = A to Z = 2 beyond a critical 
value of h. 

We first consider the case of a cylindrical slab (or disk) 
of radius R and thickness h filled with a nematic whose 
director is tangent to the two circular faces jsj (see Fig. 
IH)) . This simpler geometry captures the essential features 
of the problem and provides a suitable starting point for 
understanding thin spherical shells (see Fig. O. 




A. Slab geometry 

To estimate the energy stored in the texture of Fig. 
IHl we coarse grain the system to "blobs" of size h. The 
elastic energy arises from two sources: a long distance 
contribution from a radial texture associated with an 
5 = 1 disclination and a local energy cost for the elas- 
tic deformations inside the spherical blob in Fig. |S1 In 
the one Frank constant approximation, the former can 
be estimated as the energy nKhln (■^) of a disclination 
whose enlarged " core" of size h is given by the spherical 
blob while the latter is roughly AirKh, the energy of two 
half hedgehogs living inside the blob [SjI ■ In view of the 
azimuthal symmetry of our configuration, the director, 
n(r, z), can be parameterized by the angle 6(r, z) formed 
with respect to the S" axis of the circular slab along the 
centers of the two half hedgehogs shown in Fig. |H1 



n(r, z) — sin 9(r, z)er + cos 8(r, z)ez 



(50) 



FIG. 9: Nematic texture in a thin spherical shell. The ne- 
matic director near the pair of half hedgehogs indicated in the 
figure is well described by the one calculated for the slab in 
Fig. H 



The energy density /(G) = ^Ki (V • n)^ + \Kj, (V x n)^ 
expressed in terms of the bond angle reads 



/(e) 



Ki [smQ 



8r COS 8 — 82 sin 8 



(82 cos 8 -I- 8^ sin 8)^ 



(51) 



where K\ and ^3 are the splay and bend constants and 



8,- = ^ and 82 = ^ . In the one Frank constant 
approximation, minimization of the free energy leads to 
a non linear partial differential equation for the bond 
angle 8(r, z) 



1 a / 98 



r dr \ dr 



d^e sin 28 



dz^ 



(52) 
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The operator on the left arise from the Laplacian in cyHn- 
drical coordinates. Note there is no need to expUcitly 
consider the azimuthal angle (j) as an additional indepen- 
dent variable in view of the symmetry of the problem. 

Instead of solving this partial differential equation, we 
follow a route analogous to that in Rcf. 33], that is we 
construct an exact 2D solution for the liquid crystal prob- 
lem and then rotate it along the axis z to retrieve an 
ansatz for the 3D director configuration. The 2D solu- 
tion is the true minimum of the two dimensional liquid 
crystal elastic free energy while the 3D Ansatz does not 
minimize /(6) nor satisfy Ea. H52|l . 

To solve the 2D problem we adopt the method of con- 
formal mappings that simplifies the study of many com- 
plicated boundary problems in fluid dynamics and elec- 
trostatics 27]. Think of Fig. [Sjas a source of fluid con- 
fined to flow in a narrow and long channel {R'^ h). The 
spherical "half-hedgehog" corresponds to the source and 
the hyperbolic one at the top to the stagnation point of 
the flow. The complex potential il{w) of the desired flow 



(53) 



fl{w) — Log ( sin ( — — ) — 1 



where the complex variable is denoted by w = r + iz 
to distinguish it from the z coordinate along the axis 
of the cylinder. The velocity field is given by rt'{w) 
and the corresponding complex nematic director n{w) — 
cos B-|-i sin O is obtained by normalizing this vector field. 
Note that the prime in il'{w) denotes the derivative with 
respect to w of the complex function il = ^ + W and 
we define $7 = $ — i\l/. A straightforward but tedious 
calculation (see Appendix B) leads the functional form 
of Q{r, z) for a source at z = — -^ and a stagnation point 
at z = -|-^, namely 



tan 8(r, z) = sec ( — j sinh ( 



7rr \ 
~h) 



(54) 



This trial solution respects the boundary conditions of 
the problem, has the correct short and long distance be- 
havior in r and the expected functional form near the 
source. 

Upon inserting Q{r, z) in Ea. (|51(l and integrating the 
energy density /(O) over the cylindrical volume of the 
slab, we obtain the total elastic energy Ei stored in the 
field. The resulting energy Ei does not rely on the as- 
sumption R ^ h, as can be explicitly checked numeri- 
cally. 



El = TxKh 



In 



(55) 



where c « 4.2. Note that the functional dependence of 
El on R and h matches the expectations from the blob 
argument. Indeed the prefactor of the logarithm is uni- 
versal in the sense that it does not depend on the details 
of the trial solution near the hedgehogs, but only on its 
long distance behavior. By contrast, we expect the result 



quoted above for the coefficient c to be only an estimate 
(an upper bound) since it relies on a trial solution for 
0(r, z) that is not the absolute minimum of the free en- 
ergy. Numerical studies carried out in Ref.[31| for the 
slab geometry reported that c « 4.19. 

A competing energy minimum for the nematic is given 
by a planar texture with two s = ^ disclination lines. 
Note that -\--^ lines cannot escape in the third dimension 
|25|. A tedious but straightforward calculation shows 
that the energy Ep of the pair of disclination lines is 
given by 



--f-(-©— ^) 



(56) 



where a is a macroscopic cutoff typically of the order of 
the molecular length. These two disclination repel each 
other, and are repelled by the circular boundary, leading 
to a separation of order R. The second term of Ea. (|55|l . 
corresponding to the interaction of the two disclinations 
with the boundary and among themselves, is negligibly 
small. The third term accounts for the core energies 
of the two disclination lines 2Ec. The combination Kh 
is equal to the two dimensional coupling constant K2d- 
The relevant dimensionless ratio is -S^. 

By setting Ep = Ei we obtain the critical thickness h* 
above which the escaped "half- hedgehogs" become ener- 
getically favored compared to a single s = ^ disclination 
line: 



h* = 



(57) 



The core energy terms in Eg . 115611 reduce the previous 
estimate of c to c' = c — ^^ . A similar analysis applies 
to spherical shells which we now discuss. 



B. Extrapolation to thin spherical shells 

In principle, one could proceed along the same route 
as in the previous section and find a conformal mapping 
that provides a trial function for the bond angle <d{9,r) 
corresponding to the texture shown in Fig. |51 This is 
possible but rather cumbersome. In the case of very thin 
shells one can adapt the slab calculation by noting again 
that the energy is composed of two parts. There is a 
long distance piece arising from "combing the hair" of 
the nematic texture in the tangent plane of the sphere 
that we can read off from a suitable 2D calculation (see 
Appendix and the short distance contribution arising 
from the short distance contribution arising from the two 
pairs of half-hedgehogs at the north and south pole. 

The energy Ed of two short + 1 disclination lines placed 
at antipodal points on a sphere (the north and south pole, 
say) can be estimated by performing a 2D calculation on 
the curved surface and simply multiplying the result by 
the thickness of the layer h 



i.. = 2.i.Mln(f)-0.3+-|^^, 



-^ 



(58) 
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where the middle term accounts for the interaction be- 
tween the two dischnations in their equiUbrium positions. 
Note that this result is accurate only up to factors of the 
order of (j|) since the explicit integration over the vol- 
ume of the thin shell was bypassed. To obtain the en- 
ergy of the escaped solution, the core size a in Eq.|^ 
is rescaled to h. This will account for the integration 
of the energy density at distances of the order of a few 
h <^ R from the two hedgehogs. In these portions of the 
shell the integrand reduces to the energy density of the 
two disclination problem and hence the integration can 
be easily carried out leading the result in Eq. 158|) with a 
lower cutoff of the order h. 

The energy stored in the remaining portions of the thin 
shell is approximately given by twice the energy 4.2 T:Kh 
of the yellow blob of Fig. |S1 This estimate neglects curva- 
ture corrections of the order of (^) and arises because at 
distances of the order h the spherical shell looks locally 
like a flat circular slab as long as h ^ R. The resulting 
energy E2 of the escaped configuration reads 



E9 = 2TTKh In 



-0.3 + 4.2 



(59) 



Although the prefactor of the sub-leading term linear in 
h has only been estimated, we expect that the coefficients 
of the logarithm, which arises from large scales compared 
to h, is exact. For a spherical shell whose radius i? is a 
hundred times its thickness, the corrections from higher 
powers of -fe are indeed negligible. However for reason- 
able values of -^ , the logarithmic term of Eq. 159|) is still 
comparable in magnitude to the "sub- leading" one linear 
in h. 

The energy E4 of the tetravalent configuration can be 
evaluated using similar considerations, with the result 



Ea = TTKh ( In ( - 



0.4 



4.E, \ 
Kirh 



(60) 



Upon setting E4 — E2 we obtain the critical thickness h* 
below which the tetravalent configuration becomes ener- 
getically favored 



- e(4.i-, 



(61) 



The exponential prefactor arises from the terms linear in 
h in Eg. ljSyi) . which cannot be ignored in estimating h* 
even in the limit R ^ h. Note that an accurate deter- 
mination of the argument in the exponent would require 
knowledge of the core energies of the disclination lines. 
In fact, the exponential prefactor can be interpreted as a 
numerically significant rescaling of the core radius. 

The energy barrier between these two coexisting min- 
ima of the free energy f{Q) can be estimated by splitting 
the path connecting them in Q space in two steps. First, 
consider a continuous deformation of the escaped tex- 
ture of Fig. |51 obtained by appropriately rotating each 
nematigen until the non — escaped solution (with two 
disclinations lines of index one at the north and south 



pole) is recovered. This part of the path must be uphill 
in energy if the escaped solution was allowed to escape 
in the first place. The corresponding energy barrier AE 
is given approximately by the difference between Ed as 
calculated in Eq.ljSSJ and E2 in Eq.(|5l 



(62) 



AE = 2TrKh I In I - 1 - 4.2 
a 



The second step consists in letting each of the unsta- 
ble disclination lines split in two -fi defects and subse- 
quently separate them until they sit at the vertexes of 
a tetrahedron inscribed in the sphere. This portion of 
the path is downhill because the "non-escaped" texture 
of valence 2 is unstable. This can be proved by writ- 
ing down the energy of the pair and show that it de- 
creases monotonically as one separates them because of 
the "electrostatic- like" repulsion [i3.[l3- As a result the 
energy barrier is simply the energy difference calculated 
in Eg. lf^ . 

Upon inserting K sa 10^® dyn in Ea. H62|) and taking 
ksT w 4 lO^i'^ erg, as in Ref.jsi,!!!, we obtain 



AE _ 15h / //i\ Ec \ 

ksT nra \ \aj Knh J 



(63) 




For shells with critical thickness h* Ea. H63f) reduces to 



(64) 



where a core size of the order of 10 nm was assumed and 
the core energies were set to zero [ij . This estimate indi- 
cates that the energy barrier is very high around h* sug- 
gesting that exchange between the two minima is unlikely 
to happen by thermal activation. In a monodisperse so- 
lution of shells with thickness h, the ratio between shells 
of the two valence will be given by their Boltzman factors 
as long as equilibrium is reached. If one engineers shells 
with thickness below h* , the Z=4 configuration would be 
more likely. 



V. CONCLUSION 

We have studied the crossover from the two dimen- 
sional regime of liquid crystals confined on a spherical 
surface to the full three dimensional problem in a spher- 
ical shell. For very thin shells, the nematic ground state 
has four disclination lines sitting at the vertices of a 
tetrahedron inscribed in the ball and whose texture ap- 
proximately track the seam of a baseball. As the thick- 
ness increases, a competing three dimensional defected 
texture characterized by two pairs of half hedgehogs at 
the north and south pole becomes energetically favor- 
able. For ultra-thin shells this instability is suppressed 
and one expects a defected ground state with tetravalent 
symmetry. Estimates of the stability of this texture to 
thermal fluctuations indicate that the vibrations around 
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the equilibrium configurations of the defects should not 
be significant. The present analysis has been carried out 
primarily in the limit in which the elastic anisotropy pa- 
rameter, e = |f^<l. 

We hope to extend our investigation with a systematic 
study of the effect of elastic anisotropy on the ncmatic 
texture. It is interesting to note that in the case of pure 
bend or splay (ie. e = ±1) the ground state is given 
by only two disclinations of unit index at the north and 
south pole. This suggests the possibility that the effect 
of the elastic anisotropy may not be limited to locally 
adjusting the orientation of the director but may induce 
a change in the inter-defect interaction and hence a dis- 
tortion of the tetravalent equilibrium configuration. The 
limit of strong elastic anisotropy is also relevant to stud- 
ies of the nematic to smectic transition in a spherical ge- 
ometry for which the ratio of the bend to splay coupling 
constants -^ is expected to diverge. 

Additional experimental complications include the 
possibility of having a nematic layer of non constant 
thickness that would induce trapping of the defects in 
the regions where the layer is thinner. This effect may 
also induce a local transition to an escaped texture where 
the layer thickens in just one hemisphere so that two 
disclination lines of index i are traded for a pair of half 
hedgehogs. If that happens shells with three-fold sym- 
metry could be observed. 
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APPENDIX A: FREE ENERGY OF A VECTOR 
FIELD ON A SPHERE 

We start our analysis with the Frank free energy with 
splay and bend terms proportional to K\ and Kj, and 
expressed in terms of the covariant derivative Di n' 

F = ^Jd'^^[K\{D,nf 

+ Ks (A rij - Dj Ui) (D' n^ - D^ n')] , (Al) 



where 



D,n^ 



d,. 



ik 



(A2) 



The Christoffel connection F^^, j2j| is unchanged if the 
lower indices i and k are interchanged. As a result, the 



covariant derivative Di can be replaced by di in the sec- 
ond term of Ea. (|Al() because the covariant curl is anti- 
symmetric. It follows that 

D X n = V X n . (A3) 

The covariant form of the divergence is given by |24| 



1 



D,n' = ^d,{^i 



(A4) 



For a rigid sphere of radius R with polar coordinates 
{6*, (/)}, we have ^/g — R^sinO, F^, = — sin cos 6*, 

F ^g = F *^ = - cot e, and aU other F ^^ = 0. 

Upon adding and subtracting the expression for the 
curl of n multiplied by Ki from the first and second term 
in Ea. ljAip we obtain 

+ {K3 ^ Ki){V X n)^] . (A5) 

Similarly, upon adding and subtracting the covariant 
divergence of n multiplied by K3 from the second and 
first term in Ea. l|Al|) we obtain 

+ {Ki - K3) (V • n)2] . (A6) 

Upon adding the two equivalent expressions for F in 
Equations (|A5|) and ljA6p and dividing by two we can 
express the free energy in terms of the constants 



K 



K3-K1 
K3 + K1 
K3 + K1 



(A7) 
(A8) 



namely 



F=^ j S^^g[{D,nn{D''^ 
+e((V X n)2-(V-n)^ 



(A9) 



We now parameterize the orientation of the unit vector 
n in terms of the bond angle field a{9, 4>) that it forms 
with respect to the longitudinal direction eg which in 
polar coordinates {6, (j)) is given by 

eg = R{cos 9 sin 0, cos 9 cos (/i, — sin 9) , (AlO) 

while the orthogonal unit vector e^ is given by 

e^ = i?(— sin0sin(/>, sin6'cos0, 0) . (All) 

The components of the vector n with respect to e^ and 
eo are given by 



cos a 

IT ' 

sin a 

i? sin 61 



(A12) 
(A13) 
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Upon substituting the relevant expressions for the non 
vanishing components of the connection in the covariant 
derivative (see En. ljA2(l ') and using Ea. l)A13|l . the free en- 
ergy density f{0, 0) is given by 



4i?2 



Ki+K. 



f = ^ 



da 
'do 



e cos 2a 



2esin2a 



86 ' " 



da 

'de 



T 



where the a-dependent function Tq(0, (/>) is 
1 fda 



Tc.i9, 



smf 



COS! 



(A14) 



(A15) 



The energies of the latitudinal {a = 7r/2) and longitudi- 
nal {a = 0) tilted-molecules textures favored for e less or 
greater than zero respectively (see section Hi A|) are eas- 
ily determined by substituting the appropriate a and e 
in Ea. (|A14|l . After integrating the free energy density / 
in Ea. HA14p we obtain 



F^2nKil-\e\)l''-^de 



2ttK{1~ |e|) In 



sin(6') 



2R 



1+2^;^, (A16) 



where we have introduced a core radius, a, and corre- 
sponding core energy E^ for each defect. 

In the zero anisotropy limit (e = 0), only the first 
line in Ea. ljA15|) survives. The resulting free energy 
F — JdSf{9,(t)) then matches the one obtained using 
the spin connection in the one Frank constant approxi- 
mation upon setting Ki ~ K^ = K, 



F 



K / dSg'^{d,a - Ai){dja - Aj) , (A17) 



where dS — dddip R^ sin 9 and the metric tensor g'^'^ — 
diag {^, R'^siin^ e )' ^^^ ^^^^ "-"^ *^^ spin-connection Ai 
is the Gaussian curvature -^ [1, |23| and its only non- 
vanishing component is A^ = — cos 6. 

We now adopt a Coulomb gas representation of the liq- 
uid crystal free energy (in the one Frank constant approx- 
imation) obtained by exploiting in Ea. l|A17(l the relation 



where n(u), the defect density relative to the Gaussian 
curvature, was defined in Ea. (|A18p . The equilibrium po- 
sitions of the defects are determined only by defect-defect 
interactions because the Gaussian curvature is constant 



G = 



on an undcformed sphere. To calculate the 



Green's function r(u, u') we need to invert the covari- 
ant Laplacian defined on the sphere 



r(u,u')^- 



(A20) 



As shown below, this inversion can be accomplished by 
performing a weighted sum over eigenmodcs of the co- 
variant Laplacian, [lj|- 

We first recall that the (generalized) Green function 
r(u, u') is defined by 



Aur(u,u') = 



V9 



1 



(A21) 



where S = AwR^ denotes the area of the surface and 
A = -j=di{y^g^^dj). The presence of the second term 

on the left hand side of Ea. ljA21|) can be understood as 
follows. The Green function of the Laplacian (accord- 
ing to the usual definition without the area correction 
in Ea. (|A21|l l can be interpreted physically as the steady 
temperature response of the system to a point-like unit 
source of heat. However, for a closed system such as the 
surface of the sphere heat cannot escape. Hence, it is im- 
possible to impose a point source, that would inject heat 
at a constant rate and have the system respond with a 
time-independent distribution. To prevent energy from 
building up in such a system, we put the spherical sur- 
face of area S in contact with a reservoir that uniformly 
absorbs heat at the same rate it is pumped in. The need 
for subtracting the "neutralizing background heat" -^ in 
Ea. (|A21|l will become transparent mathematically once 
we proceed to determine r(u, u') explicitly. 

The first step consists in writing the delta function as 
a sum over spherical harmonics Y^{0, ip), 

S{u,u') _ 5{e - e')5{(j> ~ (j)') 



V~9 



R^sinO' 

^ oo m—-\-l 
1=0 m=-l 



-t'^di{d,a ~ Aj) = s(u) - G(u) = n(u) 



(A18) 



is the covariant antisymmetric tensor 
is the Gaussian curvature and s(u) = 
-j= X]i=i 9j^(u — Ui) is the disclination density with Nd 
defects of charge qi at positions u^ . The final result is an 
effective free energy whose basic degrees of freedom are 
the defect positions themselves [IJi [l3 • 



where 7*^ 
[23, G(u 



and recall the eigenvalue equation 



i?2 



/ V">' 



(A23) 



Upon substituting Ea. (|A22|l in Ea. l)A2H) and using the 
eigenvalue Equation IIA23|) , we can write down the Green 
function as 



-^UA 



jdA fdA' n(u)r(u,u')?^(u') , (A19) r(u,u')- R^^ Y. 



Yr(e,d^)Y,'' 



'; y'^yYj-'i \^ iv) ^24) 



Hl + l) 
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r!-, and used the 



We have used the fact that Yq — ,, ^^, 

neutrahzing background charge ^ in Ea. ljA21|) to cancel 
out the I = diverging mode. 

To simphfy the series in Ea. ljA24|) . we exploit the fa- 
miliar identity |35| 



771 — + / 

m— — ^ 



)Y, 



7n--^ { r\l ±J\ 



2Z + 1 

47r 



Pi(cos/3),(A25) 



where /3 is the angle (relative to the center of the 
sphere) between the directions {6*, 0} and {6'',0'} (see 
also Eq.(j3n|))- Upon substituting Ea. (|A25|l in Ea. (|A24|l . 
we find 



m— +/ 



r(u,u') 



^, G + iTT 

'I — — I ^ 



Pi{cosf3) 

47r 



(A26) 



The first term of the sum in Eq. (|A24p can be simplified 
using the following identity _36J 



■^ P;(cos/3) _ 



;=; 



In 



n 


1 + sin 




sm - 


; 


. 


\^/ J 





(A27) 



while for the second term we substitute 



E 

1=1 



PijcosP) 

l + l 



In 



In 



1 + sin 

V 


2) 


-H2J. 


-1 



with the result 

r(u,u') 



47r 



In 



1 — cos j3 



1 

47r 



(A28) 



(A29) 



Upon dropping additive constants that do not contribute 
to the energy and substituting in Ea. (|A19|l we obtain 



F 



ttK 



JVd 



Y, mj In [1 - COS A,] + E 9?^- ' (A30) 



»#i 



where the phenomenologically determined core energy 
Ec has been added by hand and reflect the microscopic 
physics not captured by our long-wavelength theory. 



APPENDIX B: LIQUID CRYSTAL TEXTURES 
AND CONFORMAL MAPPINGS 



This Appendix collects a number of results from the 
theory of complex variables relevant to the study of liquid 
crystal textures. The perspective adopted is to link the 
liquid crystal elasticity to the intrinsic geometry of the 
texture by the use of conformal transformations. The 



same method provides an elegant route to finding the flow 
lines of simple incompressible fluids in 2D |37|] and to the 
exact solution of analogous problems in electromagnetism 
and elasticity t27Ll28J. 

Nematic textures in the plane in the one Frank con- 
stant approximation can be obtained by solving Laplace 
equation, which is conformally invariant, and in complex 
coordinates {z = x + iy,^ = x — iy} reads 



dy d^a = 



(Bl) 



Here a(x, y) = a{z, z) is the bond angle that the director 
n = (cos a, sin a) forms with respect to a fixed direction, 
say the real axis x. The flat space laplacian equation 
(|B1|) is also obtained by minimizing the free energy of a 
vector field on a sphere (see Appendix El provided that 
a " stereographic projection gauge" is chosen to carry out 
the calculation 13, 38]. The stereographic projection 
maps an arbitrary point on the sphere R(^, (/>) to the 
corresponding point z — 2i?tan(|)e*'^ in the complex 
plane. The metric reads |13j 



5y 



2 1 






1 

1 



(B2) 



and the components of the gauge field in Eq. (01 are given 
by 



A^^A,^-—f^-^^ 



2iz V 1 



4_R2 



(B3) 



In this representation of liquid crystal order on a sphere, 
the Frank free energy is 



F=^ j <fz\d,a-A,W 



(B4) 



and the corresponding Euler Lagrange equation is indeed 
Ea. l|Bip . since the divergence of the gauge field is zero 
(dzAz + dzAz = 0) [lj|. The stereographic projection 
provides an example of how conformal transformations 
can be used to map physics on an arbitrary curved surface 
onto simpler planar problems. This technique can also be 
employed to analyze two dimensional order on surfaces 
oi varying Gaussian curvature (see Ref . [iTlllq 1 . 

A second use of conformal mappings is as generators 
of 2D liquid crystal textures in bounded geometries or 
in the presence of defects. Consider an analytic func- 
tion w — f{z) that maps a grid of horizontal and ver- 
tical lines in the complex plane z = x + iy onto a fam- 
ily of orthogonal curves in the w — u -\- iv plane that 
are respectively streamlines and equipotential lines of the 
corresponding flow (see Fig. I10() . Similarly, we can de- 
fine the inverse function fl(z) = f~^{z) and note that 
VL{z) = ^{x,y) + i\['(x,?/) maps equipotential lines and 
streamlines, given by the contour lines of ^{x,y) and 
^(x,j/), into a grid of vertical and horizontal lines re- 
spectively. 

The connection between liquid crystal textures and 
conformal mappings rests on the following observation: 
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complex vector field in Eq.ljBSj, with the result 



FIG. 10: (Color online) Splay rich (hyperbolic) and bend- 
rich (ellipsoidal) families of nematic "flow" lines generated by 
two s = § disclinations. The two families of flow lines are the 
equipotential and field lines of a complex function Q{z) whose 
branch cut is a horizontal line connecting the two defects. 



if the director 11(2;) forms a constant angle with respect 
to the streamlines (or the equipotential liries) of ri(z), 
then a{z) automatically satisfies Ea. ljBlll |2q |. In the 
one Frank constant approximation, the complex nematic 
director n{z) = nx{x,y) + iny{x,y) is given (up to an 
arbitrary global rotation) by 



n{z) = 



n'{z) 



(B5) 



The complex function Q'{z) denotes the derivative with 
respect to z of the complex function il — ^{x,y) + 
i"^{x, y) and we define fl = ^{x, y) — i'^{x, y). The bond 
angle is readily expressed in terms oiVt{z) via the relation 



a{z) = — 5mlog[ri'(z)] 



(B6) 



where '<sm denotes the imaginary part of a complex num- 
ber. 

As an illustration consider the simple case of two discli- 
nations on the real axis at positions x = ±1 respectively. 
The nematic director rotates by tt on a path encircling 
only one defect and by 27r on a path enclosing both (see 
Fig. I1U|) . These requirements are met by choosing the 
complex function 17 (z) = arccos(z) that is analytic ev- 
erywhere except for a branch cut on the real axis from 
a; = — ltoa; = l. The streamlines and equipotential 
lines are a family of hyperbolas and ellipses with coin- 
ciding foci at X = ±1; they correspond to two distinct 
nematic textures dominated by splay and bend respec- 
tively. The bond angle of the director oriented along the 
streamlines is easily extracted from the argument of the 



a{x,y) = arctan 



yP 



1)2 - 4a;2 



2xy 



(B7) 

In this simple example, one can explicitly check that the 
result in Ea. l|B7|l is recovered through the more familiar 
route of superposing the solutions corresponding to the 
two isolated defects 



«(.,,) 4 arctan 



V 


1 

H — arctan 
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[x- l\ 
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[x + l\ 



(B8) 



The applicability of the method of conformal map- 
pings to finding liquid crystal textures can be justified 
by means of simple geometric reasoning. We start by 
noting that the curl and divergence of a two dimensional 
vector field v, whose streamlines and orthogonal trajec- 
tories are labelled by the subscripts s and p respectively, 
can be expressed geometrically via the relations [23 



dx Vx + dy Vy 



ds \v\ + Kp \v\ 



dxVy - dyVx = V y. ^ ^ -dp \v\ + Ks \v\ 



(B9) 



(BIO) 



where Kg and Kp are the respective curvatures while ds 
and dp are the directional derivatives along the two or- 
thogonal families of level curves. The direction of in- 
creasing p is chosen to make a counterclockwise ^ angle 
with V. For example the black lines in Fig. ^| trace the 
electric field v generated by a uniformly charged plate or 
the flow lines of an ideal fluid exiting a slit of width given 
by the branch cut. Unlike the liquid crystal director in 
Eq. (jB5|) , the magnitude of v is allowed to vary with po- 
sition. By construction, such a vector field is divergence 
free and curl free, hence 



dplog\v\ 



= - (9slog|-y| 



(BU) 



(B12) 



By combining Equations (|B11|) and ljB12p . we obtain the 
geometric condition that a family of equipotential lines 
(or streamlines) needs to satisfy in order to be identified 
as level curves of an harmonic potential, namely 



dp ds 



= 



(B13) 



This condition is entirely cast in terms of the curvatures 
of the equipotential lines and streamlines without explicit 
reference to either the potential to be assigned or to the 
magnitude of the vector field v(z) [32, 113 ■ This is a 
natural language to discuss orientational order in liquid 
crystals since the director n{z) is a vector field of unit 
magnitude. 
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If we take the liquid crystal director to form a constant 
angle with respect to v(z), the curvatures k^ and Hp in 
Ea. l|B13|l can be simply cast as the directional deriva- 
tives of a along the streamlines and equipotential lines 
respectively. In fact, the curvature of these contour lines 
is the rate of change of their directions which is naturally 
parameterized by a. Hence Ea. (jB13|l reduces to 













(B14) 



The left hand side of Eq. IJB17|I is the Laplacian of a ex- 
pressed in terms of orthogonal coordinates along s and p. 
Since the Laplacian is coordinate independent, Ea. l|B14|l 
is equivalent to Ea. (|Bll) and a{x,y) represents (apart 
from an arbitrary global rotation) the desired texture 
that minimizes the Frank free energy when Ki ~ K^. 

As a byproduct. Equations IjBlip and ljB12|l can help 
to visualize how the elastic energy stored in every portion 
of the texture of Fig. ^|is distributed between bend and 
splay deformations. For most liquid crystals K^ > Ki, 
so the texture with director tangent to the streamlines 
will be energetically favored (black lines in Fig. I1U|) . In 
this case, the full Frank energy density can be rewrit- 
ten in terms of the local curvatures of streamlines and 
equipotential lines via the simple relation 



K3 (V X nf + Xi (V • n)2 = K3 Ks 



K 



1 f^p 



(B15) 



The energetically costly deformation involving bend 
takes place only around the defects in a region of radius 
of the order of their separation. Elsewhere Kg is vanish- 
ingly small. In contrast, Kp drops off more slowly at large 
distances like the inverse of the radius of a circle centered 
on the midpoint between the two disclinations. Splay de- 
formations are present throughout the system but they 
have a smaller energy cost Ki. The converse situation 
occurs if K^ < Ki so that the texture represented by the 
red equipotential lines in Fig. IIUI becomes energetically 
favored. 

The curvatures Ks and Kp are respectively the real 
and imaginary parts of the complex curvature, K{z) = 
Ks + iKp, of the mapping. This quantity can be readily 
derived from the complex potential ri(z). For calcula- 
tional purposes, it is more convenient to recast Eg. ljBlGp 



in the form k„ 



\n'\ n" 

0'2 



K{z) 



If}' If]" 



trvr) — t ■ 



n' 



(BI6) 



The reader is referred to the mathematical literature 
m nil for a proof of Ea. l|BT6|l . The intuitive signifi- 
cance of the complex curvature can be grasped by con- 
sidering how a conformal mapping f = fl~^ acts on a 
curve with local curvature k at a point in the z = x + iy 
plane. The curve is mapped onto an image curve in the 
w = u + iv plane whose curvature k' at the corresponding 
point differs from k. The curvature k' of the image curve 
is determined by two mechanisms. Firstly, the mapping 



f{z) locally stretches distances by a factor \f'{z)\, hence 
the radius of curvature of the image curve will be nat- 
urally multiplied by this amplifying factor. The second 
mechanism arises because a conformal mapping can in- 
troduce curvature even if none was originally present (in 
the isothermal net) simply by locally twisting the direc- 
tion of the isothermal net by an angle equal to arg [f'{z)]. 
For this reason, the complex derivative f'{z) is some- 
times called an amplitwist and encodes information on 
the local effect of the mapping [23 ■ The non-analytic 
function K{z) controls the amount of curvature gener- 
ated ex-novo by the mapping. For example, the mapping 
/(z) = Cos{z) transforms a grid of horizontal lines (think 
of them as a possible direction for the nematic molecules 
in the defect-free ground state) into the family of hyper- 
bolae in Fig. ^Icorresponding to a defected texture with 
two +1/2 disclinations. It is not surprising that the free 
energy density stored in the defected texture is simply 
proportional (in the one Frank constant approximation) 
to 



da 
dp 



da 
ds 



\Kiz)\' 



(B17) 



where the two elastic constants were set to be equal in 
Ea. (|B15|l and the director n(z) was parameterized in 
terms of the bond angle a{z). The Frank free energy 
is thus proportional to the complex curvature modulus- 
squared in analogy with the Helfrich free energy of a 
membrane whose derivation rests on an higher dimen- 
sional generalization of Eq. IJB15(I . 



APPENDIX C: VIBRATIONAL SPECTRUM OF 
COLLOIDAL MOLECULES 

In this appendix we provide an introduction to the 
group theoretical treatment of the vibrational spectrum 
of colloidal " molecules" . The more complicated cases of 
hexatic (p=6) and tetratic (p=4) order are analyzed in 
some detail (see Fig. [TT|l . 

The starting point of the group theoretical treatment 
of the vibrational spectrum of colloidal "molecules" is 
the observation that the defects displacements from equi- 
librium, q (see Ea. (|Hnjl . form the basis of a reducible 
representation of the point group of the molecule. If a 
molecule is acted upon by a symmetry operation, a new 
configuration will result in which the displacements of 
each defect will be permuted and transformed, but inter- 
defect distances and angles will be preserved. Here we 
take the point of view that the defects themselves are 
not permuted only their displacements, for example de- 
fect i may exchange its displacement coordinates qi with 
defect J. The liquid crystal free energy (in the one Frank 
constant approximation) is therefore invariant under all 
the operations of the point group of the colloidal defect 
array. 

The action of each operation of the group is naturally 
represented by a distinct 2N x 2N matrix {N is the num- 
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tation. Note that the characters of the original represen- 
tation are the same as the ones of the completely reduced 
one since the two differ only by a change of coordinates 
which preserves the trace. Thus, the character Xi? of the 
completely reduced representation will be the sum of the 
characters of the various irreducible representations that 
it contains 



7 



(CI) 



FIG. 11: (a) The ground state of a tetratic phase exhibits 
eight short disclination lines located at the vertices of a square 
antiprism inscribed in the sphere, (b) The twelve disclination 
that characterize hexatic order s sphere lie at the vertices of 
an icosahedron. 



TABLE II: Character for the irreducible representations of 
the icosahedral point group together with the character of 
the twenty-four-dimensional representation T generated by 
the defect displacements. 



Y 



E 



I2C5 



2OC3 



I5C2 



A 
Fi 
F2 
G 
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1 
3 
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4 
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24 
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-1 


-1 


1 








-1 


1 



2r" 



-2t 



"Note that r : 



2 • 



ber of defects) that relates the new and old defect posi- 
tions. This representation can be completely reduced by 
choosing a set of symmetry-related normal coordinates 
that are obtained from the original ones by means of 
a linear transformation. When normal coordinates are 
used, the matrixes representing the action of the sym- 
metry group can be brought in block diagonal form si- 
multaneously. Energetically degenerate linear combina- 
tion of the original coordinates form the smallest sets 
invariant under application of any symmetry operation 
of the group. The members of any one set generate an 
irreducible representation of the group. For each point 
group there is only a small number of inequivalent irre- 
ducible representations generally classified by the char- 
acters of their transformations. The characters of the 
transformations are simply defined as the traces of the 
matrices corresponding to each symmetry operation and 
they are conveniently tabulated in most texts of group 
theory [23, 113 (see Tables IIIIIII for the character tables 
relevant to the tetrahedral, icosahedral and twisted-cube 
shaped distributions of defects on a sphere). 

The task of finding the number of degenerate eigen- 
modes, n^'^\ with a given symmetry (labelled by 7) re- 
duces to counting how many times the corresponding ir- 
reducible representation appears in a reducible represen- 



where Xr labels the character of the symmetry opera- 
tion R in the irreducible representation 7. By appealing 
to the orthogonality of the characters one can write an 
expression for n^'^' in analogy with the familiar expres- 
sion for the component of a vector along a given basis 
axis Ellsi 



,(7) 



-*- Y^ (7)* 



(C2) 



where g is the number of the symmetry operations in the 
group and XR is the character of the completely reduced 
representation. Ea. (j('2|l is equivalent to Ea. (|^ quoted 
in the main text, as long as the sum over the group ele- 
ments R is replaced by a weighted sum over the classes 
in the group since the characters of group elements in the 
same class are equal. 

We now adopt the analysis of Ref.|29l l3(]| to provide 
a set of rules that produce the characters XR of the re- 
ducible representation generated by the coordinates with- 
out working out the full form of the transformation ma- 
trices. There are two key points to notice. First only 
the defects located on a symmetry axis or plane con- 
tribute to xr fl^6 trace of the transformation matrix; 
defects whose displacements are instead interchanged or 
permuted by the symmetry operation contribute only to 
the non-diagonal terms of the matrix and hence can be 
ignored in determining the character xr- Second the di- 
rections along which the displacements from equilibrium 
are measured can be chosen freely since the trace is in- 
variant upon coordinate transformations. It is generally 
convenient to choose them so that only one of the two 
displacements components is affected by the symmetry 
operation. 

As a simple example, consider a defect lying on a reflec- 
tion plane (see Fig. 112(1 . The displacement vectors before 
and after the symmetry operation is applied are {S9, S(j)} 
and {SO, —S(j)} respectively. The resulting contribution 
to the character from a single defect is thus 1 — 1 = 0. 
Inversions through a center of symmetry (coinciding with 
the center of the sphere) have a vanishing contribution to 
Xr because there are no defects there that can contribute 
to XR- 

A rotation C^ by an angle ^2^ through an n-fold axis of 
symmetry on which the defect lies leads to the following 
transformation laws for the longitudinal and latitudinal 
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TABLE III: Character for the irreducible representations of 
the tilted cube point group together with the character of the 
sixteen-dimensional representation H generated by the defect 
displacements. 
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FIG. 12: Schematic illustration of the inversion through 
an equatorial plane of symmetry (shaded circle bisecting a 
sphere) for a reflection a. The longitudinal displacement of 
the defect is unchanged while the direction of the latitudinal 
one is inverted. 



displacements 



cos 



'27^fc^ 



f2TTk'\ 



sm 




(C3) 

We measure displacements using polar coordinates with 
respect to the symmetry axis; the prime denotes the or- 
thogonal displacements after the symmetry operation C^ 
is applied. Inspection of Ea. ljC3|) shows that the contri- 
bution from C^ to the character xr is equal to 2 cos (^^) 
times the number of defects lying on the axis of rotation. 
On the other hand the contribution to XR from the im- 
proper rotation S^ is zero. To see this note that the 
symmetry operation S^ is a rotary reflection achieved by 
performing a successive rotation through an (alternating) 
axis followed by a reflection in the plane perpendicular 
to the axis k times. An example of a molecule possessing 
the symmetry operation 5*4 is methane, CH4, with the 
carbon atom lying at the intersection between an alter- 
nating axis and the reflection plane (See Fig. 5-2 in Ref. 
^3|)- The tetrahedral defect configuration considered 
in this paper does not posses any defect at the position 
occupied by the carbon atom of the methane molecule. 
More generally, the possibility of having a defect whose 
equilibrium position is unchanged by the rotary reflec- 
tion is ruled out because such defect would have to lie 
off the spherical surface at the intersection between the 
alternating axis and the plane of reflection. 

To sum up, each of the characters, xRj of the com- 
pletely reduced representation formed by the displace- 
ment coordinates is given by the number of atoms whose 
equilibrium positions are not changed by the symmetry 
operation R times its fundamental character as derived 
in the previous paragraphs. Similar results that apply to 
unconstrained molecules whose atoms have three dimen- 



sional displacements are listed in Table 6-1 of Ref. |29l |. 
The resulting characters for the tetrahedral, icosahedral 
and tilted cube defects configurations are listed in Ta- 
bles IIIIIll Upon using Ea. (|43() and the character table UTI 
we can decompose the 24 dimensional representation, Y, 
formed by the displacements from an icosahedral equilib- 
rium configuration into irreducible representations. The 
result reads 



Y = 2H + 2G + 2Fi 



(C4) 



The three rigid body rotations correspond to one of the 
two triplets in Fi while the remaining 21 independent 
normal coordinates form energetically degenerate multi- 
plets with the following degeneracy factors: 2 quintets, 2 
quartets and 1 triplet. This analysis is confirmed upon 
direct diagonalization of the representation Y which leads 
the 21 non- vanishing eigenvalues Ai, with the multiplici- 
ties shown bold in parenthesis 



A, = {0.87X (5), 0.09X (5), 
0.74X (4), 0.22X (4), 0.96x (3)} , 



(C5) 



Note that the normal modes in the second quintet are 
much "softer" than the rest. 

A similar analysis applied to the twisted cube configu- 
ration of defects leads to the following decomposition of 
the (defects' displacements) representation, S, 



E^2Ei + 2E2 + 2E^ + A1+A2 + B1+B2 



(C6) 



where the three rigid body rotations are contained in 
one of the two doublets E3 and in the singlet A2, which 
leaves five doublets and three singlets for the eigenvalues 
(i. Direct diagonahzation of S leads four doublets, one 
triplet and two singlets of non- vanishing eigenvalues, 

Q = {1.37X (3), l.aix (2), 0.89X (2), 

0.47X (2), 0.06X (2), 0.11, 1.26} (C7) 

The discrepancy between the degeneracies found by di- 
rect diagonalization on one hand and group theory on 
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the other is caused by an accidental symmetry of the po- 
tential energy of the tilted-cube arrangement of defects. 
Hence the first triplet is to be interpreted as the missing 



doublet and singlet that happen to have the same energy 
even if there is no symmetry reasons to expect so. The 
modes in the last doublet of Ea. (|C7|l are the softest. 
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